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The quantum antiferromagnetic spin- 1/2 Ising model on a triangular lattice and analogous fully 
frustrated Ising model on a square lattice with quantum fluctuations induced by the application of 
the transverse magnetic field are studied at finite temperatures by constructing an exact mapping 
onto a purely classical model with a more complex interaction. It is shown that in weak fields the 
temperatures of the phase transitions separating the critical phase from the ordered and disordered 
phases in both models are proportional to the magnitude of the field. 
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I. INTRODUCTION 

In recent years the interest to quantum frustrated mag- 
netic systems is constantly growing due to the perma- 
nent appearance of new materials belonging to this class 
(see Ref. 1 for a recent review). In many cases, an inter- 
play between fluctuations (quantum and thermal) and 
a macroscopic degeneracy of the classical ground states 
makes understanding the properties of such systems a 
nontrivial task. 

The best known example of a frustrated magnetic sys- 
tem is the antiferromagnetic Ising model on a triangular 
lattice, in which each triangular plaquette has to contain 
at least one frustrated bond with higher than minimal en- 
ergy. It is well known from the exact solution^ that this 
classical model is disordered at any positive temperature 
T > 0, whereas at T = it is characterized by a finite 
residual entropy and an algebraic decay of correlations.'^ 
A quantum analog of this model can be defined by the 
Hamiltonian, 

ur) J 

where J > is the coupling constant of the nearest- 
neighbor interaction, T > is proportional to the magni- 
tude of the transverse magnetic field inducing quantum 
fluctuations, and (t^ are the Pauli matrices, and the 
summation in the first term is performed over all pairs of 
nearest neighbors {jj') on a triangular lattice. The classi- 
cal antiferromagnetic Ising model is recovered in the case 
r = 0. 

Besides being applicable for the description of an easy- 
axis antiferromagnet with spin 1/2 and a triangular lat- 
tice, the quantum antiferromagnetic Ising model defined 
by Eq. (1) is also of interest as a representative of a 
wider class of the fully frustrated transverse-field Ising 
models, in which coupling constants can have both signs 
(Jjji = ± J), but on each plaquette of the lattice the num- 
ber of the antiferromagnetic bonds (with Jjj' = J > 0) is 
odd.^ The quantum fully frustrated Ising models on var- 
ious lattices have been investigated^"^° mostly in view 
of their relation to the quantum dimer models^^ with 
vanishing potential energy term. A method for con- 



structing more general models implementing a contin- 
uous crossover between a frustrated Ising model and a 
dimer model (on the dual lattice) whose Hamiltonian in- 
cludes both the kinetic and potential terms has been pro- 
posed recently in Ref. 12. 

The first investigation of the finite-temperature phase 
diagram of the quantum antiferromagnetic Ising model 
on a triangular lattice was undertaken by Isakov and 
Moessner.-'^'^ Relying on the analogy with the classical 
antiferromagnetic Ising model on a layered triangular 
lattice"'^^ in a system with a finite size in the direction per- 
pendicular to the layers they argued that at < F < Fc 
(where Fc oc J is the critical value of F at zero tem- 
perature) this model has to have two finite-temperature 
phase transitions of the Berezinskii-Kosterlitz-Thouless 
type with a critical phase between them (like in the clas- 
sical six-state clock model) and confirmed these conclu- 
sions with the help of numerical simulations. Soon after 
that, Jiang and Emig^^ proposed a derivation (based on a 
renormalization group analysis) demonstrating that the 
temperatures of these two transitions and T2 should 
behave themselves as 

T{=,2CxFln''(Fc/F), (2) 

where ly is the critical exponent of the three-dimensional 
XY model and T^/T[ = 9/4. 

However, the prediction of Ref. 15 was based on the 
renormalization-group analysis assuming the possibility 
of describing the evolution of the system in Euclidean 
time in the framework of the continuous approximation. 
This approach requires the size of the system in the time 
direction /? = 1/T to be much larger than the typi- 
cal time between spin flips (inversely proportional to F), 
which at temperatures given by Eq. (2) is fulfilled only 
for Fc - F < Fc. 

The main aim of the present work is to study the be- 
havior of Tf and in the quantum antiferromagnetic 
Ising model on a triangular lattice at F ^ J, where the 
continuous approximation is not applicable. Our result 
consists in finding that in this range of parameters both 
transition temperatures are just proportional to F and 
do not diverge with the increase of J as predicted by 
Eq. (2). We also demonstrate that in the fully frustrated 
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transverse-field Ising model on a square lattice the situa- 
tion is the basically the same: the ordered and disordered 
phases are separated by the critical phase with both tran- 
sition temperatures being proportional to the magnitude 
of the field in weak fields. Up to now, the investigation 
of this model has been focused on its zero-temperature 
properties. 

II. ANTIFERROMAGNETIC MODEL ON A 
TRIANGULAR LATTICE 

A. The case of the infinite coupling constant 

We start by considering the quantum antiferromag- 
netic Ising model on a triangular lattice in the special 
case of the infinite coupling constant, J = oo. In such a 
case, the Hilbert space of the model is restricted to the 
states with only one frustrated bond in each triangular 
plaquette (that is, the ground states of the classical anti- 
ferromagnetic Ising model labelled below by index a) and 
their linear combinations. The finite-temperature parti- 
tion function can be then written as 

a 

with 

Wa = {a\exp{-pH)\a) (4) 

and /3 = l/T. 

In the infinite-temperature limit (/3 0) one gets 
Wa — > 1, that is, Z is reduced to the zero-temperature 
partition function of the classical antiferromagnetic Ising 
model. At /3 > weights Wa become dependent on the 
structure of the state a. A convenient way for describ- 
ing this dependence consists in introducing an effective 
classical Hamiltonian Hcff(fl) defined by the relation 

Wa = cM-'Hofiia)] . (5) 

At /?r <C 1 the spin flips induced by the second term in 
Eq. (1) are rare which allows one to calculate the sta- 
tistical weights Wa corresponding to different states per- 
turbatively. In such a way, one obtains that in the lowest 
order in (3T, Hesia) is just proportional to M(a), the 
number of spins in state a each of which can be flipped 
without taking the system out of its Hilbert space, 

n,s{a) = -^iPT)Hlia) . (6) 

The terms of the fourth and higher orders in /3 depend on 
the structure of the state a in a more complex way, but 
at /3r <C 1 the main role is played by the terms included 
into Eq. (6). 

It is well known that the set of the ground states of the 
classical antiferromagnetic Ising model on a triangular 
lattice can be mapped onto the states of a solid-on-solid 



(SOS) model describing [111] facet of a cubic crystal. -'^^ In 
terms of the SOS representation, each Ising spin aj — ±1 
is replaced by integer height variable hj in such a way 
that the relation 

h,-h,^-l±^^[tl (7) 

is satisfied for every pair of neighboring sites on the trian- 
gular lattice. Eq. (7) assumes that one of the three basic 
vectors of the triangular lattice (whose sum is equal to 
zero) is directed from site j to site j' and not from j' to j 
(which would correspond to the opposite sign of hji ~hj). 

When each triangular plaquette contains only one frus- 
trated bond (with OjUji = -1-1), Eqs. (7) unambiguously 
define the values of all integer variables hj as soon as the 
value of hj is chosen for one of the sites. On the other 
hand, any state of the SOS model with hji — hj = +1, —2 
corresponds to some ground state of the antiferromag- 
netic Ising model. 

It follows from the known properties of the classical 
antiferromagnetic Ising model on a triangular lattice that 
when all allowed states of such an SOS model enter the 
partition function with the same weight, this model is in 
the rough phase, in which correlations of heights diverge 
logarithmically,^'' 

K 

9jk = {{hj - hkf) w — \ogrjk , (8) 

where rjk is distance between sites j and k and 
K = Ko = 18. 

The same surface representation can be used for inter- 
preting the quantum spin model (1) with J = cxd as a 
quantum SOS model in which the amplitude of height 
jumps, hj — >■ /ij ± 3 is given by T. Alternatively, the use 
of Eq. (4) allows one to speak about the classical SOS 
model in which now the weights corresponding to dif- 
ferent configurations of the surface are no longer equal 
to each other and depend on the configuration. Since in 
terms of the height representation flippable spins corre- 
spond to points where the surface has no local slope, in 
terms of the SOS representation functional M{a) is a 
measure of the flatness of the surface. Then it is clear 
from the negative sign in Eq. (6), that the interaction 
described by this equation suppresses the fluctuations of 
the surface. 

It has to be emphasized that the classical SOS model 
defined by the Hamiltonian 'Hoff(a) is exactly equivalent 
to the original quantum spin model. The price one has to 
pay for the reduction of a quantum model at a finite tem- 
perature to a purely classical one is that it is impossible 
to write down the explicit form of T-Lcs{a), Eq. (6) being 
applicable only at /3T <^ 1. However, a number of prop- 
erties of the classical SOS model defined by Hoff (a) can 
be discussed without knowing the exact form of 'Heff(a)- 

In particular, it is clear that the decrease of T (increase 
of /3 = 1 /T) leads to the suppression of the factor K in 
the correlation function (8) and finally has to induce a 
phase transition into the ordered (flat) state. Note that 
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at T = 0, a surface described by a quantum SOS model 
always has to be in the ordered (flat) phase^® with a well- 
defined value of (h) and saturation of correlation function 
gji; at large distances. 

However, for /3T <C 1 the suppression of factor K has 
to be small, whereas the phase transition to the ordered 
state will take place when this factor is suppressed from 
Kq = 18 down to Kf = 4 (see Ref. 17). Note that this has 
to be so independently of whether in the ordered phase 
the spins on all three sublattices are magnetically ordered 
or if on one of them the average magnetization is zero 
(see discussion in Ref. 13 and Ref. 15). In terms of the 
SOS representation the first case corresponds to having 
(h) integer and the second one to having (h) half-integer, 
but in both cases the periodicity of the ^.-dependent ef- 
fective potential is the same, and the critical value of K 
is determined by this periodicity.^^ 

Apparently such a pronounced suppression of K can- 
not happen while /3r <^ 1 and requires T ^ T, whereas 
at T ^ r the system has to remain in the critical phase. 
One can conclude that at J = oo the phase transition 
between the critical and the ordered states of the orig- 
inal spin model takes place at a finite temperature 
proportional to F whose value follows from the relation 
K (Tf ) = Kf. 



B. The case J < oo 

When J < oo, the Hilbert space of the model is sub- 
stantially extended because now all configurations of the 
Ising spins aj = ±1 are allowed. In terms of the SOS 
representation, this corresponds to the appearance of the 
possibility of the creation of screw dislocations^'' on go- 
ing around which the integer variable hj interpreted as 
height instead of returning to the same value changes by 
so-called Burgers number b = ±6. Each dislocation is cen- 
tered around a plaquette containing not one but three^*^ 
frustrated bonds with Cjaji = -1-1 which therefore can be 
identified with the dislocation core. In the framework of 
the path integral description of a quantum system, dislo- 
cations are linear topological excitations and have either 
to form closed loops in space-time or to cross the whole 
system in the direction of the Eucledean time. 

In terms of the original spin variables the closed dis- 
location loops in space-time correspond to the spin fiips 
which are prohibited at J = cx). At F <C J such pro- 
cesses can be taken into account in the framework of the 
perturbation theory. The most important of them is the 
second-order process which leads to the decrease of the 
systems energy by the amount proportional to F^/J per 
each spin which cannot be flipped without increasing the 
energy of the system. At T ~ F ^ J this gives a correc- 
tion to K of the order of F/ J and therefore leads only to 
a small shift of Tf with respect to its value at J = oo. 

On the other hand, the dislocations crossing the whole 
system in the time direction can lead to the disordering 
of the critical phase. At F/J — 0, when the system does 



not experience evolution in the Eucledean time these dis- 
locations can be identified with that of the classical SOS 
model. The core energy of such dislocations is propor- 
tional to J, which makes their fugacities at T ^ J ex- 
ponentially small. The main difference which appears at 
F/ J <^ 1 is that quantum fluctuations lead to a small 
negative correction to the core energy but this is irrele- 
vant for further reasoning. 

When logarithmic interaction of dislocations is strong 
enough, they are bound in neutral pairs and the sys- 
tem has the same properties as in the absence of 
dislocations.^^ The strength of this interaction is deter- 
mined by the same parameter K as the amplitude of the 
fluctuations of h. The dislocations with Burgers number 
b remain bound in neutral pairs only for 

K<K^, = j (9) 

whereas at if > K2 there appear free dislocations whose 
proliferation leads to the disordering of the critical phase. 
In the antiferromagnetic Ising model on a triangular lat- 
tice dislocations have b = ±6 and, accordingly, A'l = 9 
(as found in Ref. 17). 

Since K2 = 9 is smaller than Kq = 18 but larger than 
Kf = 4, one can conclude that at F <C J, the disso- 
ciation of dislocation pairs leading to disordering of the 
critical phase takes place at temperature T2 ^ F which is 
higher than Tf and therefore there exists a finite interval 
of temperatures < T < where the system remains 
critical. The value of is determined by the relation 
K{T^) = = 9. As Tf , the transition temperature T| 
is basically proportional to F and remains finite when J is 
taken to infinity at a finite F with the difference between 
the values of at J = 00 and at F <C J < 00 being ex- 
ponentially small in J/F. Note that there are no reasons 
to expect the ratio T^/Tf to be equal to K2IK1 = 9/4 
(as it was proposed in Ref. 15) because factor K is not 
proportional to T but depends on the ratio F/T in a more 
complicated way. 



III. FULLY FRUSTRATED MODEL ON A 
SQUARE LATTICE 

The same approach can be applied to the fully frus- 
trated transverse-field Ising model on a square lattice. In 
the case J = 00, the Hilbert space of any fully frustrated 
transverse-field Ising model is defined by the set of the 
ground states of the classical fully frustrated Ising model 
on the same lattice, which is isomorphic to the full set of 
states of the classical dimcr model on the dual lattice."* 
On a square lattice, the classical dimer model without 
any interaction of dimers is exactly solvable^"'^ and allows 
for a mapping onto a SOS modeP^ with integer height 
variables hj defined on the sites of the same square lat- 
tice. For a given configuration of Ising spins Cj, the values 
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of hj can be defined by the relation analogous to Eq. (7), 
hf -h,^-{l + 2T,,,a,a,, ) = { ^3 (10) 

where Tjj> = Jjj'/J = ±1- Eq. (10) assumes that the 
square lattice is divided into two equivalent square sub- 
lattices (A and B) and site j belongs to sublattice A and 
site j' to sublattice B. In the opposite case the sign of 
hji — hj would be the opposite. 

When each square plaquette contains only one frus- 
trated bond (with Jjj'Crjajr = J > 0), Eqs. (10) unam- 
biguously define the values of all integer variables hj as 
soon as the value of hj is chosen for one of the sites. 
On the other hand, any state of the SOS model on a 
square lattice with hji — hj = +1,-3 corresponds to 
some ground state of the fully frustrated Ising model on 
the same lattice. 

Since the classical system of nonintcracting dimers on a 
square lattice is in the critical phase, ^'^ the corresponding 
SOS model is in the rough phase with the asymptotic 
behavior of the height-height correlation function given 
by Eq. (8), where according to Ref. 22 K = Kq = 32. 
Like in the case of a triangular lattice, the description 
in terms of the classical SOS model in which all allowed 
configurations enter the partition with the same weight 
is applicable to the fully frustrated transverse-field Ising 
model with J = cx) in the classical limit (3 — (that is, 
T = 00). 

Having T < 00 decreases the fluctuations of the sur- 
face. In the lowest order in /3T the Hamiltonian of the 
classical SOS model describing the system is again given 
by Eq. (6). In terms of the dimer representation, this 
Hamiltonian corresponds to having an attraction be- 
tween parallel dimers belonging to the same plaquette. 
The classical dimer model with such an interaction has 
been studied in Ref. 24. In our system, the effective 
Hamiltonian has this form only at the highest temper- 
atures, whereas with the decrease of temperature the in- 
teraction of the more distant dimers also starts to play 
a role. However, since variables hj are integer, the phase 
transition to the ordered phase takes place at the tem- 
perature at which factor K is decreased down to = 4 
(exactly like in the case of a triangular lattice), inde- 
pendently of what is the exact form of the dimer-dimer 
interaction. In accordance with that, at J — 00 the tem- 
perature of the phase transition to the ordered phase is 
proportional to F and for F ^ J < 00 the correction to 



the value of this quantity is small. 

The second phase transition is related to the appear- 
ance of free (unpaired) dislocations crossing the whole 
system in the direction of Euclidean time. Like in the 
case of a triangular lattice, dislocations can be associ- 
ated with the plaquettes which contain not one but three 
frustrated bonds. For K < K2 [where i^f is given by 
Eq. (9)] they are bound in neutral pairs which dissoci- 
ate when the value of K reaches K!^- In the SOS model 
describing the ground states of the fully frustrated Ising 
model on a square lattice Burgers numbers b arc equal 
to ±8 and, accordingly, A'2 — 16 (see Ref. 24). Since 
K2 = 16 is smaller than Kq = 32 but larger than = 4, 
one can again make a conclusion that the critical phase 
exists in a finite interval of temperatures < T < T2, 
whereas the phase transition leading to the disordering of 
the critical phase takes place at temperature T2 at which 
K{T^) = ATI = 16. For F < J this temperature like Tf 
is proportional to F. 



IV. CONCLUSION 

In the present work wc have demonstrated that in the 
quantum antifcrromagnctic transverse-field Ising model 
on a triangular lattice, as well as in the fully frustrated 
Ising model on a square lattice, the temperatures of the 
phase transitions separating the critical phase from the 
ordered and disordered phases in weak fields are propor- 
tional to the magnitude of the field. The analysis of Jiang 
and Emig^^ leading to a different conclusion fails in weak 
fields (F <C J) because the continuous approximation 
used in Ref. 15 for the description of the spin fluctua- 
tions in Eucledean time requires the size of the system in 
the time direction /? to be much larger than the typical 
time between spin flips which is inversely proportional 
to F. For T ^ T[2, this condition is fulfilled only when 
Fc — F <C Fc, whereas out of this range the continuous 
approach cannot be trusted. The results of the numerical 
simulations of the antifcrromagnctic model on a triangu- 
lar lattice^'^ are consistent with Tf and T2 being propor- 
tional to F at small F. 

The author is grateful to F. Mila and S. Wenzel for a 
useful discussion. 
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